Study of gossamer superconductivity and antiferromagnetism in the t-J-U model 
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The d-wave superconductivity (dSC) and antiferromagnetism are analytically studied in a renormal- 
ized mean field theory for a two dimensional t-J model plus an on-site repulsive Hubbard interaction 
U. The purpose of introducing the U term is to partially impose the no double occupancy constraint 
by employing the Gutzwiller approximation. The phase diagrams as functions of doping 8 and U 
are studied. Using the standard value of t/J = 3.0 and in the large U limit, we show that the 
antiferromagnetic (AF) order emerges and coexists with the dSC in the underdoped region below 
the doping 8 ~ 0.1. The dSC order parameter increases from zero as the doping increases and 
reaches a maximum near the optimal doping 8 ~ 0.15. In the small U limit, only the dSC order 
survives while the AF order disappears. As U increased to a critical value, the AF order shows up 
and coexists with the dSC in the underdoped regime. At half filing, the system is in the dSC state 
for small U and becomes an AF insulator for large U. Within the present mean field approach, We 
show that the ground state energy of the coexistent state is always lower than that of the pure dSC 
state. 

PACS: 74.25. Jb, 71.10.Fd, 74.72.-h, 74.25. Ha 
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I. INTRODUCTION 

In spite of tremendous theoretical and experimental 
efforts dedicated to the studies of the anomalous prop- 
erties of high T c superconductors (HTS), a full under- 
standing of these materials is still far from the final stage. 
As a basic point, it is known that much of the physics 
should come from the competition between the d-wave 
superconductivity (dSC) and antiferromagnetism. Ex- 
perimentally, it is generally suggested that the ground 
state evolves from the antiferromagnetic (AF) state to 
that of the dSC order as the carrier density increases 1 . 
However, since the early days of HTS, there also have 
been persistent reports of the coexistence of the dSC and 
AF orders 2-8 in various cuprate samples. Especially in 
the recent neutron scattering experiments, the commen- 
surate AF order has been observed in the underdoped 
superconducting YBa2Cu3C>6.5, providing the unambigu- 
ous evidence for an unusual spin density wave state co- 
existing with superconductivity (dSC) 6 . Therefore it is 
necessary to develop a microscopic theory in which both 
the antiferromagnetism and the dSC are treated equally 
in order to understand the ground state property of the 
cuprate superconductors. 

Theoretically, it has been widely accepted that the es- 
sential physics of cuprates can be effectively described 
by the two dimensional Hubbard model or its equiva- 
lent t-J model in the large U limit 9,10 . Using the varia- 
tional Monte Carlo (VMC) method, several groups pro- 
posed wave functions with coexisting AF and dSC orders 
and found that the coexisting state has a lower energy 
than either the pure dSC order or the pure AF state in 



the underdoped regime 11-14 . Although the slave particle 
mean field theory for the t-J model was originally in- 
troduced to investigate the formation of the RVB state 
or the superconducting order 9 ' 15-18 , it also has been ap- 
plied to study the coexistence of the dSC and AF orders 
in this system 19 ' 20 . Stimulated by the idea of the "gos- 
samer superconductors" proposed by Laughlin 21 , Zhang 
and co-workers 22 employed the t-J-U model with the 
Gutzwiller projected wave function 23 to investigate the 
superconducting order parameter and the electron pair- 
ing gap (or the RVB order parameter). There 22 the on- 
site Coulomb interaction U is introduced to partially im- 
pose the no double occupancy constraint for the strongly 
correlated electron systems. In the large U limit, their 
result 22 is consistent with that of Kotliar and Liu 16 using 
the slave boson mean field approach for the t-J model. 

Following Ref. [22], we report a further investigation 
of the same model by taking the AF order explicitly into 
consideration. Within the Gutzwiller renormalized mean 
field theory, we find that for large Coulomb repulsion 
U, there is a coexistence between AF and dSC orders 
below the doping level S ~ 0.1. The coexisting state 
always has a lower energy than that of the pure dSC 
state. The dSC order parameter increases from zero as 
the doping increases in the underdoped region and then 
reaches a maximum near the optimal doping 5 ~ 0.15, 
after that it decreases to zero at S ~ 0.35 with increasing 
doping. When the magnitude of U is reduced, the AF 
order parameter decreases very quickly with increasing 
doping, and the coexistent region is squeezed toward low 
doping regime until it disappears for U < 5.3t, where 
the "gossamer superconductivity" is found even at half 
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filling. 

The paper is organized as follows. In Sec. II, we out- 
line the theoretical framework. The t-J-U model is in- 
troduced and the Gutzwiller variational approach is for- 
mulated. A renormalizcd Hamiltonian is obtained and 
further studied within the mean field theory. In Sec. Ill, 
our numerical results are displayed and compared with 
those from other theories and experiments. In Sec. IV, 
a summary of the paper will be given. 



II. THEORETICAL FRAMEWORK 

We start from the t-J-U model on a square lattice 22 , 
H = H t + H s + H v , (1) 

with 

H t = -t^2(ClC i+na + h.c), 

ifja 
ifj 

Hu = U^h t] h ti , (2) 

where r) = x and y, C\ a {Cia) is the electron creation 
(annihilation) operator, S, = ^Z aa , C^BaaiCia' /2 is the 
spin operator with a = (a X7 <j y ,<j z ) as Pauli matrices, 
hia = Cj a Cicr, U is the on-site Coulomb repulsion, t is 
the hopping parameter, and J is the exchange coupling 
constant. In the Hamiltonian (1), the U term is intro- 
duced to partially impose the no double occupancy con- 
straint. In the limit U — > oo, the model is reduced to the 
t-J model. 

To study the Hamiltonian (1) with the Gutzwiller vari- 
ational approach, we take the trial wave function as 

|V>> =P G |Vo(A d ,A Q/ , M )), (3) 

where Pq is the Gutzwiller projection operator and it is 
defined as 



P G = n i [l-(l- ff )n iT J 



■ill 



(4) 



here g is a variational parameter which takes the value 
between and 1. The choice g — corresponds to 
the situation with no doubly occupied sites(£7 — > oo), 
while g = 1 corresponds to the uncorrelated state(£/ = 
0). Vo(Ad, A a /, /j,)) is a Hartree-Fock type wave func- 
tion, where A<f,,A /,/j are the parameters representing 
dSC, antiferromagnetism and chemical potential, respec- 
tively. The nature of |?/>o} depends on the expected long 
range behavior. Since it is the purpose of this paper 
to study the interplay between antiferromagnetism and 
dSC, we will adopt the wave function which includes both 
the dSC and antiferromagnetism in a unique variational 
space 12 ' 13 . 



With help of the trial wave function (3), the variational 
energy E var = (H) is given by 

E var = ^}, H \! P) = NUd+{H t ) + (H s ), (5) 



where 



(V>|V) 

{Ht) = 
(H s ) = 



<V>|V> : 



(6) 



N is the total number of the lattice sites and d = (ri^n^) 
is the average double occupation number. Obviously, the 
double occupancy can be modulated by U. 

In the calculation of the variational energy we adopt 
the Gutzwiller projection method which was formulated 
originally for the Hubbard Hamiltonian. A clear and sim- 
ple explanation 24 was given by Ogawa et al. and by Voll- 
hardt. In their scheme, the spatial correlations are ne- 
glected, and the effect of the projection operator is taken 
into account by the classical statistical weight factors. In 
this way the hopping average and the spin-spin corre- 
lation in the state | ijj) are related to those in the state 
| tpo) through the following relations 

■4 



g a {4>o \ Si • Sj |Vo). (7) 



(1>\Si- S, | V) 



In the thermodynamic limit, one has 24 

2 _ d{l-n + d) (n-2wr) 2 
9 ~ (1 - r){l - w)wr (n-2d) 2 ' 



(8) 



and the renormalization factors can be derived to have 
the following expressions, 



9t = 



n-2d 
n — 2rw 



1 — w)(l — n + d) lw ~ 
1 — r V r 



(l-r)(l-n + d) , /_r rf 



1 — w 



, n — 2d . o 

( ) • 

n — 2wr 



(9) 



(10) 



Here n is the average electron number per site. In order 
to consider the AF order, the square lattice is divided 
into two sublattices A and B. For sublatticc A we assume 
(hii) = r = t£ +m and (n^) = w = ^—m, i.e., a net mag- 
netization +to at each site. For sublattice B the electron 
occupation numbers r and w are exchanged, meaning 
the magnetization —m at each site. Here m represents 
the AF order parameter in the state \il>o). These renor- 
malization factors, g t and g s , quantitatively describe the 
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correlation effect of the on-site repulsion. We will further 
comment on this point below. 

In terms of these renormalization factors, the varia- 
tional energy E var — (H) is rewritten as 

E var = (H ef/ ) , (11) 
where H e ff is the Gutzwiller renormalizcd Hamiltonian: 
Heff = 9tH t + g s H s + Hjj 

= -g t tY J {Cld+^ + h.c.) 



+ gsjJ^Si ■ S i+fj + NUd. 



(12) 



if} 



In the mean field approximation the renormalizcd 
Hamiltonian (12) can be rewritten as 



H MF = NUd + -Ng s J(A 2 + X 2 ) + 2Ng s Jm 2 



-aA af (Cl a C k+Qa + Cl +Qa C ka )} 
- A d rj k {C-k[Ck^ - C- k+ QiC k+ Qi 



(13) 



where the electron chemical potential /i has been added, 
Q = (tt, it) is the commensurate nesting vector, and the 
prime on the summation symbol indicates that k is lim- 
ited to half of the Brillouin zone. In the above equation, 
we have introduced respectively the electron pairing or- 
der parameter, the hopping average and the staggered 
magnetization 



i+ril 



— Ci-tC, 



= A (—A) when rj = x (y) , 
m=(-l) i (C^C n -C} l C il ) /2 , 



(14) 
(15) 
(16) 



with -f k = 2(cosk x + cosky), rj k = 2(cosk x — cosk y ), 
e k = -{gtt+lg s Jx)lk, A d = §c/ s JA, and A af = 2g s Jm. 
Here the parameter A d is always associated with the fac- 
tor r\ k in Eq. (13), which implies that the supercon- 
ductivity has a d-wave like symmetry. The mean field 
Hamiltonian (13) is easily diagonalizcd, giving rise to four 
bands, ±E\ k and ±E 2k with 



Ei k = yfijTk - A 4 ) 2 + i&dVk) 2 , 
E 2k = V(-^-li) 2 + (A dr]k ) 2 , 

^ = ^t + Kf ■ (17) 

Here A^r/h and A a f can be regarded respectively as the 
energy gap associated with the dSC and the AF order 
parameter. The ground state energy is given by 



E var /N = Ud- n& - + E 2k ) 



+ -g s J(A 2 + X 2 ) + 2g s Jm 2 . (18) 

By minimizing the ground state energy, we can ob- 
tain the self-consistent equations for the quantities A 
(the electron pairing order parameter), \-> m (staggered 
magnetization), d and the chemical potential it as follows 



(19) 



X 



m 



^ 7fe a ( "E 1 T + ^- ) ' (20) 



2 k 



1 WA„/ A-ju -Cfe-M> 



2N ^ Cfe E lk 



E 2k 



1 ^ dEvar dg t | dE var dg s ^ 



= UN 



ANg s J dg t dm dg s dm 
dE var dgt dE var dg s 



dg t dd dg s dd 



1 v^'/^fe - A* , -jfc ~ M 
E 2k 



3;< 



+ 



■)• 



(21) 
(22) 
(23) 



For each doping S, all the parameters A, \, m, d and 
/j, are determined self-consistently by the Eqs. (19)-(23). 



III. RESULTS AND DISCUSSION 

Now we summarize our results. Firstly we discuss the 
average double occupation number d as a function of U. 
Our calculated results at the doping 6 = 0.0 (solid line), 
0.05 (dashed line) and 0.1 (dotted line) for the parameter 
t/J = 3.0 at the temperature T=0 are shown in Fig. 1. 
We find that the average double occupation number d 
at S = 0.0 decreases linearly as function of U till U = 
9.M, where d shows the similar behavior of discontinuity 
as reported in Ref. [22]. But for the doped cases, our 
numerically obtained d as functions of U do not show this 
discontinuity, and they become flattened and decrease 
slowly at large U. 

The Gutzwiller renormalization factors gt and g s as 
functions of doping S for the parameters t/J = 3.0 and 
U = 20t at T = are shown in Fig. 2. The dashed 
lines are the corresponding results when the AF order is 
not considered or m is fixed to zero. As we mentioned in 
Sec. II, these factors quantitatively reflect the partially 
enforced no double occupancy constraint due to the on- 
site Coulomb repulsion U. For large U, the effect of the 
Gutzwiller projector operators is to reduce the kinetic en- 
ergy and enhance the spin-spin correlation. We find that 
at low doping, the AF order suppresses the magnitude of 
g s while g t is only slightly affected. 
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FIG. 1. The average double occupation number d as a func- 
tion of U at doping S — 0.0 (solid line), 0.05 (dashed line), 
and 0.1 (dotted line) for the parameter t/J — 3.0 at T=0. 
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FIG. 2. The Gutzwiller renormalization factors gt and g s 
as functions of doping 5 for the parameters t/J — 3.0 and 
U = 20t at T — (solid lines). The dashed lines are the 
corresponding results when the AF order is not considered, 
i.e., m is fixed to zero. 

In Fig. 3, we plot the self-consistently obtained order 
parameters A and m as functions of doping 5 for the pa- 
rameters t/J = 3 and U = 20t at T = 0. The dashed 
line is the corresponding A when the staggered magne- 
tization m is set to zero. It should be noticed that these 
parameters are the expectation values under the wave 
function IV'o)- It is clear that the electron pairing order 
parameter A is drastically suppressed at low doping by 



the AF order. At half filling, A is reduced to zero and 
to reaches to its maximum value. Near <5 ~ 0.1, the AF 
order vanishes while A shows a peak. 




FIG. 3. The self-consistent parameters A and m as func- 
tions of doping 5 for the parameters t/J — 3.0, U = 20t at 
T = 0. The dashed line is the A with m is set to zero. 

We now discuss the dSC order parameter Asc and AF 
order parameter tuaf under the wave function \ip) 1 which 
are defined as 

= A S c (-Asc) when r\ = x (y) , (24) 
ttiaf = (-l) 1 ^^ - C*Cii>/2 • (25) 

In the Gutzwiller approximation, these parameters are 
easily obtained from A and m with the following renor- 
malization factors: 

Asc = .9aA, 

m AF = g m m. (26) 
Similar to the method of deriving g t and g s , we obtain 

2 



9A 



n~2d 
2(n - 2rw) 



(l-w){l-n + d) 



+ 



+ 



(l-r)(l-n + d) 



9m = 



1 — w 



n-2d 
n — 2wr 



+ 



(27) 



(28) 



In Fig. 4 we plot the dSC order parameter Asc> AF 
order parameter itiaf and the electron pairing gap (or 
the RVB order parameter 22 ) A^ = |g s JA as functions 
of doping S for t/J = 3.0 and U = 20t at T = 0. From 
this phase diagram, we find that the AF and dSC or- 
der parameters coexist for a wide doping range, up to 
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S ~ 0.1, in the ground state. It can also be seen that 
the AF order parameter is a monotonically decreasing 
function of S, but the dSC order parameter shows a non- 
monotonic dome shape: it increases from zero as the 
doping increases in the underdoped region and then has 
a maximum near 8 ~ 0.15, after which it decreases to 
zero at 8 ~ 0.35 with increasing doping. Although the 
present approach applies only at T = 0, the supercon- 
ducting transition temperature T c (8) is expected to ex- 
hibit a similar 8 dependence, and to have a maximum at 
the optimal doping 8 ~ 0.15. It should be noticed that 
the electron pairing gap Ad is also reduced to zero at 
half filling because of the presence of the AF order. This 
is quite different from the case in Rcf. [22], where the 
AF order is not considered, and the electron pairing gap 
increases as the doping decreases. 




FIG. 4. The dSC order parameter Asc, AF order param- 
eter rtiAF and the electron pairing gap Ad as functions of 
doping 6 for U = 20t and t/J = 3.0 at T = 0. 

In order to further understand the effect of the 
Coulomb repulsion U on the ground state behavior, cal- 
culations for several other values of U are performed. In 
Fig. 5, we plot the calculated results for U = 5t, 7t, lOt 
and 15t with t/J — 3 and T = 0. It is clearly seen that 
with decreasing U, the AF order decreases very quickly 
with increasing doping, and the coexistent region of the 
AF and dSC orders is squeezed toward lower doping. Par- 
ticularly for U = 5i, the coexistence disappears, and the 
AF order is completely suppressed by the prevailing dSC 
order. To illustrate more clearly the dependence of the 
order parameters on U , we present the parameters Asc 
and rriAF as functions of the Coulomb repulsion U for 
doping 8 = 0.0(a), 8 = 0.05(b) and 8 = 0.1(c) at T = in 
Fig. 6. At half filling (see Fig. 6(a)), for small Coulomb 
repulsion U < 5.3t, only the dSC order persists. As U 
increases up to U — 5.3t, the AF order begins to show 
up and coexists with the dSC and the transition appears 
to be a second order. At U = 7t, there is a discontinuity 
in the slope of vhaf and the dSC order gets completely 
suppressed by the AF order at U > It where our sys- 



tem becomes an AF insulator. For U > 9.3i, the double 
occupancy number d drops discontinuously to zero. As 
a result, the magnitude of rriAF jumps from 2.7 to 3.8 
and becomes U independent for large U. And with in- 
creasing doping (see Fig. 6(b)), the AF order exists only 
for larger U while the dSC order is always in presence. 
But for doping 8 > 0.1 (see Fig. 6(c)), the AF order 
completely disappears independent of the magnitude of 
U. 
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FIG. 5. The dSC order parameter Asc, AF order param- 
eter m,AF and the electron pairing gap A d as functions of 
doping 8 for different values of U with t/J = 3.0 and T = 0. 
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FIG. 6. The dSC and AF order parameters Age and vtiaf 
as functions of the Coulomb repulsion U for different dopings 
with t/J = 3.0 and T = 0. 
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With the help of these self-consistent parameters, let us 
compare the ground state energy obtained from Eq.(18) 
with that of Ref. [22] in which the contribution from the 
AF order was neglected. In Fig. 7, we plot our ground 
state energy E var /N as a function of doping 8 using the 
parameter t/J = 3.0 for several different values of U 
(see the solid lines) . The dashed lines here correspond to 
the results when the contribution from the AF order is 
not included, i.e., m is fixed to zero 22 . From Fig. 7, we 
conclude that the ground state energy with the AF order 
considered is always lower than that without it. 
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FIG. 7. Doping dependence of the ground state energy for 
several different U for the parameter t/J — 3.0. The dashed 
lines are the corresponding results when the AF order is not 
considered, i.e., m is fixed to zero. 

We now discuss the relevance of our calculations to 
other theories. Although the t-J model, derived from 
the large U Hubbard model, was originally introduced 
to study the superconductivity based on the RVB theory 
without AF order 9 ' 15-18 , the inclusion of the AF order 
based on the same approach was done at a much later 
stage. In all these studies, the no double occupancy con- 
straint has been globally enforced. Using the t-J or a 
similar model and based upon other type of mean field 
approximations, there exist several works 19 ' 20 ' 25-28 in- 
vestigating the existence of both AF and dSC orders in 
the system. While the double occupancy is globally ex- 
cluded from the standard t-J model, our current t-J-U 
model with finite U allows partial enforcement of the no 
double occupancy constraint, and to understand the sub- 
tle effect due to the electron-electron correlation. For the 
case of small U, our results show that only the dSC order 
exists in the ground state, which describes the physics of 
the " gossamer superconductor" . In the limit of infinite 
U, the t-J-U model is reduced to the t-J model. In this 



case our phase diagrams show that the AF and dSC or- 
ders coexist with each other from small 5 up to 5 ~ 0.1, 
and after that the AF order completely disappears. This 
feature is in good agreement with the VMC results for the 
t-J model 12-14 . At the same time, we notice that the co- 
existence between the AF and dSC orders persists up to 
optimal doping S ~ 0.15 in the slave-boson scheme 19 ' 20 . 
We would mention that the similar large coexistence can 
be obtained if we neglect the derivatives of gt and g s with 
m in our derivation of the self-consistent equations, i.e., 
replace Eq. (21) with the following one, 



m 



1 wA a/ g fc 
2N^ & 1 



-6c - /" 



E 



ik 



2 k 



(29) 



In this way, we can perform similar calculations as above. 
In Fig. 8, we present such a phase diagram with t/J — 
3.0 and U — 15t at T = 0. It can be seen that in this 
case, the AF and dSC orders coexist up to doping 5 ~ 
0.18. But it seems that such a large coexistent region is 
not favored by the experimental and simulation results. 
Moreover, based on this approximation, the system at 
half filing would always be an AF insulator, independent 
of the magnitude of U. This is contrary to what has been 
obtained from our current approach based on minimizing 
the total energy of our system. 
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FIG. 8. The dSC order parameter Asc, AF order param- 
eter vtiaf and the electron pairing gap as functions of 
doping 8 for U = 15t and t/J = 3.0 at T = 0. Here the 
derivatives of g t and g s with m in the self-consistent equa- 
tions are neglected. 

So far the experimental evidences for the coexistence 
of the AF and dSC orders in cuprate superconduc- 
tors seem not conclusive. For example, the long range 
AF order observed in the insulating L^-^Sr^CuC^ is 
sensitive to doping 1 , which disappears rapidly at x ~ 
0.03. But there also existed several experimental re- 
sults which appeared to indicate the coexistence of anti- 
ferromagnctism and superconductivity over a wide dop- 
ing range in cuprate superconductors 2-8 . Especially, 
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the AF order was claimed to have been observed in 
underdoped YBa2Cu3C>6.5 and YBa2Cii306.6 supercon- 
ductors by neutron scattering experiments from differ- 
ent groups 6 ' 8 , ft is apparent that more experiments are 
needed to confirm the coexistence of the long range AF 
order with the dSC state in HTS. 



IV. SUMMARY 

In summary, we have studied the coexistence of the an- 
tiferromagnetism and dSC in a renormalized mean field 
theory based on the Gutzwiller approximation for a two 
dimensional t-J-U model. The role of the Hubbard inter- 
action U is to partially enforce the no double occupancy 
constraint, and it provides us with a better understand- 
ing of the subtle effect due to the electron-electron corre- 
lation. Our results show that the AF and dSC orders co- 
exist below the doping 5 ~ 0.1 at large U with t/J = 3.0. 
And we find that the coexisting state has a lower ground 
state energy than that of a pure dSC state. The dSC 
order increases from zero as doping increases in the un- 
derdoped regime and reaches a maximum near the opti- 
mal doping S ~ 0.15, after which it decreases to zero at 
S ~ 0.35 with increasing doping. With decreasing U, the 
coexistent region is squeezed toward low doping. There 
is no coexistence between AF and dSC orders for small 
U(< 5.3i), where the AF order is completely suppressed 
and the "gossamer superconductivity" is found even at 
half filling. For the large U, our system at half filling is 
always an AF insulator in which both the electron pairing 
gap and the dSC order parameter are suppressed to zero. 
Our result at large U should correspond to the physical 
regime. The reason why the existence of the long range 
AF order has not been firmly confirmed by experiments 
in the underdoped HTS is probably due to the neglecting 
of the AF fluctuations in the mean field approximation. 
It is believed that the effect of the AF fluctuations may 
break the long range AF order into short range orders, 
and this conjecture needs to be examined more carefully 
in future theories and experiments on cuprate supercon- 
ductors. 
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